%Compute the probability of a AND b
function p = prob_joint(a, b, m, n, alpha)
pa = prob_one(a, m, alpha);
pb = prob_one(b, n, alpha);
p = pa * pb;
end